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ABSTRACT 

Wc study dynamical interactions of star-planet binaries with other single stars. We derive analytical 
cross sections for all possible outcomes, and confirm them with numerical scattering experiments. We 
find that a wide mass ratio in the binary introduces a region in parameter space that is inaccessible to 
comparable-mass systems, in which the nature of the dynamical interaction is fundamentally different 
from what has traditionally been considered in the literature on binary scattering. We study the 
properties of the planetary systems that result from the scattering interactions for all regions of pa- 
rameter space, paying particular attention to the location of the "hard-soft" boundary. The structure 
of the parameter space turns out to be significantly richer than a simple statement of the location 
of the "hard-soft" boundary would imply. We consider the implications of our findings, calculating 
characteristic lifetimes for planetary systems in dense stellar environments, and applying the results 
to previous analytical studies, as well as past and future observations. Recognizing that the system 
PSR B1620-26 in the globular cluster M4 lies in the "new" region of parameter space, we perform a 
detailed analysis quantifying the likelihood of different scenarios in forming the system we see today. 
Subject headings: stellar dynamics — planetary systems — globular clusters: general — methods: 
ri-body simulations — pulsars: individual (1620-26) 



1. INTRODUCTION 

The subject of planets in dense stellar systems such 
as globular clusters came to life with the discovery that 
the binary pulsar PSR B1620-26 in the globular clus- 
ter M4 is orbited by a distant companion of approxi 



ter M4 IS orbited by a distant companion oi approxi- 
mately Jupiter mass BackeHlT993l iThorsett et al .11199^ 



iSigurdsson et alJ2003l andreferences therein) . The pres- 
ence of such a system near the core of M4 is, at first, 
puzzling for several reasons. For example, the metallic- 
ity of M4 is low (O.OSZq), making the formation of such 
a planet via the standard core accretion model very un- 
likely, an d sparking the cr eation of alternative formation 
theories (jBeer et al.l2004. and references therein). In ad- 
dition, the planet is on such a wide orbit about its host 
that the system should have been destroyed by dynam- 
ical interactions on a very short timescale (~ 10^ yr) in 
the moderately high density environment of M4's core, 
implying that it must have been formed recently via an 
alternative mechanism such as a dynamical interaction 
involving a binary. 

Campaigns searching for "hot Jupiter" planets (with 
orbital periods < 10 d) in 47 Tuc have yielded null detec- 
tions, instead placing upper limi ts on the hot Jupiter 
frequ ency there (jCilliland et alj 120001 iWeldrake et al.l 
l2005j) . It now appears that the lack of hot Jupiters 
detected in 47 Tuc is probably more closely connected 
with the low metallicity of the environment, and not 
simply due to destruction via dynamical encounters in 
the high density stellar environment. This result is con- 
sistent with the strong correlation between host star 
metallicity and planet frequency in observed systems 
ijFischer fc Valentill2005j) . Searches for planets via tran- 
sits in the higher metallicity environments of open clus- 



ters are currently underway^. 

A key ingredient in the formation and evolution of 
planetary systems in dense stellar environments that is 
absent in lower density environments is dynamics. If ei- 
ther the stellar density is large enough or the planetary 
system is wide enough, it is likely to undergo a strong 
dynamical scattering interaction with another star or bi- 
nary, perturbing its properties or possibly destroying it in 
the process. The importance of dynamics for wide mass 
ratio binaries cannot be overstated, especially when one 
considers the wealth of low-mass objects such as brown 
dwarfs that are known to exist in open and globular clus- 
ters, and may be in binaries. In addition, the importance 
of dynamics as early as the stage of star formation is 
now becoming clear. There is now convincing evidence 
that a distant encounter with another star is responsi- 
ble for the currently observed prop erties of distant So- 
lar System objects such as Sedna llKenvon fc BromlevI 
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|2^ iMOTbidelii fc Levisonl2004t iKobavashi et allfoO 

The object HD 188753 is a hierarchical triple star system 
composed of a main-sequence star in an « 12.3 AU orbit 
about a compact main-sequence binary. What makes 
this system interesting is that the single star hosts a 
hot Jupiter in orbit abou t it, with a semimajor axis of 
w 0.04 AU ()Konackil|200,5t) . The current configuration of 
the system is puzzling, since the binary companion would 
truncate the initial disk out of which the planet could 
have formed at a radius much smaller than that required 
by either the core accretion or gravitational fragmenta- 
tion paradigms. It is now clear that the current configu- 
ration of HD 188753 must have resulted from a dynami- 
cal scattering interaction i n a moderately dense stellar 
envir onment fPfahl .20051 iPorteeies Zwart fc McMillad 
I2ti05l) . 

The types of dynamical interactions we have been 
* see |http://star-www.st-and.ac.uk/~kdhl/transits/table.html| 
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discussing have been studied in great detail and vol- 
ume since at le ast the niid-1970s for the case of com- 
parab le masses (|Hillsl 119751 iHeggidllQTSt iHeggie fc 
120031 and references therein). The case of disparate 
masses in the binary (such as for a planet ary system) 
has been treated in less detail l|Hilla.l984: .Hills fc Disslvl 
Il989t iSigurdssonI Il992() , but with considerable atten- 
tion paid to the survivability of planetary systems in 
particular iLaughhn fc Adam£lll998HBonnell et"aI]l200H 
iDavies fc Sigurdssonll200l|K To the best of our knowl- 
edge, there has not yet been a comprehensive treatment 
of dynamical scattering interactions with wide mass ra- 
tio systems in the literature, as previous work has typi- 
cally been concerned with a numerical treatment of a cer- 
tain aspect of the problem. This paper is an attempt to 
remedy this situation by carefully treating a wide range 
of parameter space both analytically and numerically, 
and identifying physically distinct regions of the space, 
with an emphasis on application to physically interesting 
problems. 

In particular, the questions wc would like to answer 
in this paper are the following: for a dynamical inter- 
action between a star-planet system and another star, 
under what conditions does the interaction result in: 1) 
preservation of the planetary system, in the sense that 
the planet remains in orbit about its parent star, even if 
the orbit is perturbed; 2) exchange of the incoming star 
for the planetary system's host star; 3) destruction of 
the planetary system, in the sense that the planet is not 
bound to either star after the interaction; and 4) physical 
collisions involving two or more of the three objects? For 
the outcome in question 2, we would also like to know 
the resulting properties of the planetary system's orbit. 
In particular, we would like to know the regimes in which 
the planetary system orbit shrinks and those in which it 
widens. 

In the language of binary scattering, the interaction of 
a binary planetary system with a single star can be writ- 
ten symbolically as '[Mi m] M2', where the Mj represent 
the stars, and m represents the planet. The possible out- 
comes can be enumerated succinctly as shown in Tabled 
The outcome "pres" is preservation of the planetary sys- 
tem, and can result either from a strong interaction or a 
weak fly-by. The planetary system can have shrunk or 
widened in the process. The outcome "ex_s" is an ex- 
change in which the host star of the planetary system is 
exchanged for the incoming star. The outcome "ex_p" 
is an exchange in which the two stars remain bound, 
ejecting the planet. The outcome "ion" leaves all ob- 
jects (stars and planet) unbound from each other. The 
remaining outcomes listed in the table are various per- 
mutations involving one or more collisions. The outcome 
of a stable hierarchical triple is omitted from the t able 
since it is classically forbidden (|Heggie fc Hutil200^ . It 
is possible to be formed as a result of a binary-single in- 
teraction if tidal effects are present. However, we ignore 
such effects in this paper. 

Clearly, case 1 above corresponds to outcome "pres", 
case 2 corresponds to outcome "ex_s", and case 3 cor- 
responds to outcomes "ex_p" and "ion". The outcomes 
corresponding to case 4 are the remaining listed in the 
table. In this paper we address the questions posed above 
by calculating cross sections for the processes listed, as 
well as characteristic lifetimes. In section |21 we treat as 
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Fig. 1. — Plot of the critical velocity Vc and orbital speed Horb 
for a planetary system composed of a IMq star and a 10~^Mq 
planet in an orbit with semimajor axis a, encountering a IMq 
star with relative velocity at infinity Vao- Dotted lines represent 
characteristic velocity dispersions for open and globular clusters, 
and orbital sizes for hot Jupiters and Jupiters at 5AU, as shown. 



much of the problem as we can analytically. In section 
|31 we confirm our analytical predictions with numerical 
experiments. In section0]we determine under what con- 
ditions the planetary system will shrink or expand as the 
result of an exchange encounter, and the average amount 
by which it shrinks or expands. In section [S] we discuss 
the implications of the results in the context of current 
literature and observations. In section El we quantify the 
likelihood of the possible dynamical formation scenarios 
of the planetary system PSR B1620-26 in the globular 
cluster M4. Finally, in sectional we summarize and con- 
clude. 

2. ANALYTICAL CONSIDERATIONS 

As described above, we symbolically label the inter- 
action between a planetary system and a single star as 
'[All fn] -M2'. We also use the symbols Mj and m for 
the masses of the bodies, along with a for the planetary 
system semimajor axis, e for its eccentricity, foo for the 
relative velocity at infinity between the planetary system 
and single star, and b for the impact parameter. If one 
treats all bodies as point particles, and only considers 
their masses for the time being, there arc only two veloc- 
ity scales in the problem. The first is what has histori- 
cally been called the critical velocity, which is the value 
of Voo for which the total energy of the binary-single sys- 
tem is zero, and is given by 



\ MO 



1/2 



(1) 



where /i = (Mi -I- rnjK'h/ {Mi + M2 + m) is the reduced 
mass of the binary-single star system. The second is the 
orbital speed of the planet relative to its parent star, and 
is given by 



I'orb 



G(Mi + m) 



(2) 



for a circular orbit. Looking at cqs. ^ and l(2Jl, wc see 

that if all the masses in the problem are comparable, then 
i^c ~ Vorh- If, on the other hand, one of the members of 
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TABLE 1 

Possible outcomes of a dynamical scattering interaction between a 
planetary system, '[mi m]', and a single star, 'm2'. 



outcome 


description 


nickname 




Ml m] M2 
M2 m] Ml 
Ml M2] m 
Ml M2 m 


preservation 
exchange host star for incoming star 
exchange planet for incoming star 
ionization 


"pres" 
"ex_s" 
"ex_p" 
"ion" 


[Mi-.m Mj\ 
Mi-.m Mj 
[Mi-.Mj m] 
Mi-.Mj m 


colhsion of planet with one star, resulting in binary 
collision of planet with one star, resulting in ionization 
collision of two stars, resulting in binary 
collision of two stars, resulting in ionization 




Mi-.Mf.m 


collision of all objects 





Note. — The symbol denotes a physical collision. 



the binary system is much less massive than the other 
member (as is the case for a planetary system), and we 
take Ml = M2 = M, and q = m/M ^ 1. then we can ap- 
proximately write Vc ~ ^/qvorh- For a typical planetary 
system, q ^ 10"'^, and the two velocities differ by more 
than an order of magnitude. Fig. Q] plots the two veloci- 
ties for a planetary system with q = 10~^ and M = Mq 
as a function of a. In essence, introducing a mass ratio in 
the binary opens up a strip in parameter space (between 
the two lines in Fig. ^ in which the physical charac- 
ter of dynamical scattering interactions is fundamentally 
different from that of the comparable-mass case. 

In the literature on binary scattering involving stars of 
comparable mass, Vc is the velocity delin iiting the bor - 
der between "hard" and "soft" binaries l)Heggielll975|) . 
If Voo < Vc-, the binary is in the hard regime, and will 
on average harden (a will get smaller) as a result of an 
encounter. On the other hand, if Wqo > the binary is 
in the soft regime, and will on average soften (a will get 
larger) as a result of an encounter — that i s, if the binary 
is not dest royed comp letely. However, iHills fc Disslvl 
()1989(1 and iHillsl ijlQQOj) suggest that when considering 
binaries with disparate masses, the boundary is more ac- 
curately given by Voo = Uorb, and not by Voo = Vc, and 
suggest the use of the terminology "fast-slow boundary" 
rather than "hard-soft boundary" , since it is the rela- 
tive speed of the binary members that is physically rel- 
evant. In the literature considering the survivability of 
planetary systems (and high mass ratio systems in gen- 
eral) in dense stellar environments, both Vc and Vorh have 
been used as the single characteristic velocity delimiting 
the boundary between ionization on the high-velocity 
side and hardening on the low- velocity side. However, 
as we will see below, for planetary systems the hard- 
soft botmdary lies at Vao = Vc, while the characteristic 
timcscalc for the survivability of a planetary system in a 
dense environment drops markedly only at Voo = Worb- 

Having written down the two characteristic velocities 
in the problem, thus dividing parameter space into three 
distinct regions (as shown in Figure the question 
naturally arises: What is the physical character of the 
scattering interactions in each of the three regions? For 
Voo < Vc, the total energy of the binary-single system 
is negative. We thus expect that for any "strong" en- 
counter (defined approximately as one in which the clas- 
sical pcricenter distance, rp, between the single star and 
the planetary system is of order a), the interaction will be 



resona nt, in the sense that i t will survive for many orbital 
times (jHeggie fc Hutll2003|) . It will behave as a small 
star cluster, most likely ejecting its lightest member — in 
this case the planet — yielding the outcome "cx_p (res)". 
For Tp > a, the likely outcome is a "direct" exchange 
of the incoming star for the planet ("ex_p"). We thus 
expect (Tcx.p = o-ox.p(rcs) + o-ox.p(non-ros) to bc thc dom- 
inant cross section for Voo < Vc- For Voo > Vovh, the 
interaction is impulsive, since the timescale of the inter- 
action is a small fraction of the binary orbital period. We 
thus expect that any "strong" interaction will ionize the 
system, making aion the dominant cross section in this 
regime. For v^, < Vqo < forb the timescale of the interac- 
tion is comparable to the orbital period, so the outcome 
is not a priori obvious, and must bc treated numerically. 

Below, as far as it is possible, we calculate using ana- 
lytical techniques the cross sections for each outcome in 
each of the three regions of parameter space. 

2.1. Voo < Vc 

First, since for < Vc the total energy of thc binary- 
single system is negative, it is clear that ionization is 
classically forbidden. Thus 

CTion = . (3) 

We should also mention that for the outcome of preser- 
vation the only method we have for determining if an 
interaction was "strong" is whether or not it was reso- 
nant. Thus it should be clear that for every region of 
parameter space in Figure ^ 

^prcs(non— res) ^ ; (4) 

since every very distant passage of the single star by the 
b inary preserves it vi a a non-resonant interaction. 

iHeggie et all ()1996|) considered interactions of "hard" 
binaries (in the sense that Voo < Vc) with single stars 
for a wide range of possible mass ratios in the binary 
and between the binary members and thc incoming single 
star. Using analytical techniques, they calculated the 
scaling of the cross sections for each possible outcome. 
They then fit to numerically calculated cross sections to 
determine the weakly mass-dependent "coefficients" on 
each cross section. From their eqs. (13) and (14) with 
m2 ~ m3 (since Mi w M2 here), we sec that 

'^cx_p(rcs) ~ '^cx_p(non— res) ■ (5) 
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q 

Fig. 2. — Plot of the ratio of the cross section for resonant 
exchange in which the planet's original host star is ejected, to that 
in which the planet is ejected, o"ex_s(res)/''"ex_p(rea) i s-s a function of 
the binary mass ratio q = m/M. Here we have assumed that both 
stars have mass M. Note how, as q is decreased from unity, the 
test-particle limit is reached relatively quickly, at q 0.1. 



The cross section for non-resonant exchange of the planet 
for the incoming star is given by their eq. (17) (divided 

by 2 since 0'ox_p(rcs) ~ ''^cx_p(non — res))- 




where Mt ~ il/i+m+A/2 is the total mass of the binary- 
single system, the exponential term is the "coefficient" fit 
to the numerical results, and the amn are given in their 
Table 3. 

The cross section for direct exchange of the incoming 
star for the host star is also given by their eq. (17), with 
the appropriate permutation of masses: 




The resonant cross section for this process, crex_s(rcs)i can 
be calculated approximately by statistical techniques. 
Since for a resonant interaction we expect the three-body 
system to behave like a small star cluster, we can as- 
sume that the less massive body (the planet) and the 
more massive bodies (the stars) reach energy equiparti- 
tion, yielding the approximate equality 

ww^,™ = Mvl^j^i , (8) 



where Vc,m is the one-dimensional velocity "dispersion" 
of the planet, v^.m is the one-dimensional velocity disper- 
sion of the stars, and we let M\ = M2 = M . Assuming 
that the three body system reaches thermal equilibrium 
on a timescale shorter than the interaction timcscale, 
each species obeys a Maxwellian distribution: 

/(.,n,..) = 4..^^^-|^e-''V-S, (9) 

where n is the number density of the species, Va its one- 
dimensional velocity dispersion, and v its speed. To find 
the ratio of the probability for ejecting one of the stars to 
ejecting the planet, one simply integrates the velocity dis- 
tribution for each species from the escape velocity, Dcsc, 
to infinity and considers the ratio. Using the fact that 
Vcsc ~ Sv^Wctj where « v^^m is the one-dimensional 
velocity dispersion of the full three-object system, we 
have 

O- ex.p (res) S^^^ f {v , Tl, Va / y/q)dv 

26-6^1 +erfc(^/6) 
2e-69y^ + erfc(y6g) ' 

where erfc(z) = 1 — erf (2) is the complement of the er- 
ror function. Note that there is no factor of two in the 
numerator since the two stars in the system are distin- 
guishable, in the sense that we distinguish between "ex_s 
(res)" and "pres (res)". This function is plotted in Fig. El 
Note how, as q is decreased from unity, the test parti- 
cle limit is reached relatively qui ckly, at q 0.1. Th e 
same qualitative result is found in lFregeau et alJ l)2004() . 
as shown in their Figure 10. 

Since in a resonant interaction the system loses all 
memory of its initial configuration, we expect that for 
Ml = AI2 = M, the probability of preserving the system 
(and ejecting the incoming star) in a resonant interaction 
is equal to that of ejecting the host star and keeping the 
incoming star: 

^prcs(rcs) '-^cx_s(rcs) ■ (H) 

Finall y, we remark on the li kelihood of collisions in this 
regime. iFregeau et al.l l)2004l) have studied in detail the 
probability of collisions during binary interactions, in- 
cluding binaries with very small mass ratios, and found 
that for systems in which the radius of the low-mass 
member is independent of its mass (as is approximately 
the case for planets), the test particle limit is reached at 
q fa 0.1 (see Figure 10 of that paper). In other words, 
the collision cross section is independent of the mass of 
the low-mass member below q w 0.1. Furthermore, the 
collision cross section is roughly as large as the strong 
interaction cross section. Therefore, we expect the com- 
bined cross section for any outcome involving a collision 
to be on the same order as crox.p(ros)- This estimate says 
nothing about the relative frequency of the different out- 
comes involving collisions, which requires numerical cal- 
culations to infer, and which we do not treat in this pa- 
per. However, we do consider collisions briefly in the 
discussion of PSR B1620-26 in sectional 

2.2. Voo > Worb 
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In the regime Voo > forb, since > Vc, the total 
energy of the system is positive, so resonant scattering is 
forbidden: 



^prcs(rcs) 



cx_s(rcs) 



^cx_p(rcs) 



0. 



(12) 



Since Woo > I'orb, the timcscale of the scattering inter- 
action is typically much smaller than the orbital period 
(unless rp ^ a). In this case the interaction is in the im- 
pulsive limit, and ca n be handl ed quite efficiently with 
analyt ical techniques l)Hutlll983|) . From eq. (5.1') of lHuti 
(119831). the ionization cross section is 



40 ^ Mi 1 

TTflG — 

3 Ml + m 



(13) 



Exchange in the impulsive regime occurs when the in- 
coming star, M2, passes much closer to one member of 
the binary than the other, and undergoes a near-180° 
encounter in which all of its momentum is transferred to 
that member. In order for M2 to remain bound in the 
binary, it must receive nearly zero recoil in the frame of 
the remaining binary member, requiring that the mass 
of M2 be nearl y equal to the mass of the binary mem- 
ber it displaces (|Hutlll98'^ . Thus exchange of the planet 
for the incoming star, crex_p(non-rcs)j is vanishingly small 
unless M2 ~ m. Exchange of the planet's host star for 
the incoming star, on the other hand, is more likely since 
the stars' m asses are m ore likely to be comparable. From 
eq. (5.2') of lHutI l)1983j) . the cross section for exchange of 
the planet's host star for the incoming star is given by 



cx_s(non— res) 



—n-G^M\M + m)- 

3 a V, 



(14) 



with Ml = M2 = M. 



2.3. Vc < Voc < Vorh 

In the regime Vc < I'oc < '^orb, the interaction is nei- 
ther resonant nor impulsive, so the techniques used above 
cannot be used here. For q 1 the interaction is essen- 
tially the hyperbolic restricted three-body problem, but 
since the timescale of the interaction is of order a few 
orbital periods, the standard techniques of perturbation 
theory are of limited use. Instead, one must appeal to 
numerical methods to discern the nature of the interac- 
tions. It is possible, however, to make order of magni- 
tude estimates of the relevant cross sections. Whenever 
M2 passes within Vp < a oi the binary, its force on the 
planet is comparable to that of Mi, in which case the 
orbit of the planet should be sufficiently perturbed so as 
to unbind it from the host star and either eject it from 
the system or leave it bound to Al2. Thus we estimate 
that the total cross section CTion + o"ox^(non-ios) should be 
roughly equal to the strong interaction cross section. In 
the absence of any more detailed information, we simply 



assume a\. 



-^cx_s(noii — res) ; 



yielding 



GM, 



^cx_s(non— res) 



ira- 



(15) 



since gravitational focusing dominates. Since in 
this regime the planet is effectively a test particle, 
o'ox.p(non-ros) IS vanishingly small. 






pres (res) 





ex_s (res) 


A 


ex_s (non-res) 


D 


ex_p (res) 




ex_p (non-res) 


A 


ion 




Fig. 3. — The dimensionless cross sections, a/na^, for the out- 
comes involving no coUisions in Table ^ as a function of Voo / Vc 
for scattering interactions between a planetary system with mass 
ratio q = m/M = 10~^ and eccentricity e = 0, and a star of mass 
M. Points plotted are as described in the legend. Circles repre- 
sent resonant interactions, while triangles represent non-resonances 
( "direct" interactions) . The error bars plotted are one-sigma. The 
vertical dotted lines are the critical velocity and orbital speed. The 
solid lines are the analytically predicted cross sections from section 
1^ while the dashed line is the order of magnitude estimated cross 
section from the same section. 



3. NUMERICAL CONFIRMATION 

To confirm the analytical predictions of the previous 
section, we have calculated the cross sections by numer- 
ically integrating many scattering interactions. We have 
used the Fewbody numerical toolkit to p erform the scat- 
tering interactions ijFreeeau et al.ll2b04|) . For calculating 
the cross sections, we use a general i zation of the tech- 
nique described in lMcMillan fc HuO ljl996j) in which for 
each annulus in impact parameter we perform only one 
scattering interaction, instead of the many required by 
their technique. This allows us to more precisely calcu- 
late the uncertainties in the cross sections, in particu- 
lar yielding accurate asymmetric error bars. T o handle 
the small number statistics, we use the tables in lGehrel^ 
©86). 

FigureOlshows the dimensionless cross sections, a/ira'^, 
for the outcomes involving no collisions in Tabled as a 
function of Voo/vc for a planetary system with mass ra- 
tio q = 10-3 and e = 0, and with Mi = M2 = M. 
Points plotted are as described in the legend. Circles 
represent resonant interactions, while triangles represent 
non-resonances ("direct" interactions). The error bars 
plotted arc one-sigma. The vertical dotted lines are the 
critical velocity and orbital speed. The solid lines are the 
analytically predicted cross sections from section|2 while 
the dashed line is the order of magnitude estimated cross 
section from the same section. We have calculated the 
same cross sections using St arlab ijMcMillan fc Hutll996t 
iPortegies Zwart et al.ll200if) . and the agreement is excel- 
lent. Note that with the exception of sectional all calcu- 
lations presented in this paper are performed with e = 0. 
Previous results show that most cross sec tions of interest 
do no t depend strongly on eccentricity l)Hut fc BahcaU 
Il983j) . Indeed, we have calculated the cross sections in 
Figure |31 using a thermal eccentricity distribution and 
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find only a slight increase in the resonant cross sections 
(and corresponding decrease in the non-resonant cross 
sections) relative to the e = case. 

Looking first at the region Voo > f orb , we sec dion plot- 
ted in blue and (Tcx_s(non-rcs) plotted in red, along with 
the corresponding analytically predicted cross sections 
(cqs. ^21 ^nd respectively) plotted as solid lines. 
The agreement is clearly quite excellent. 

Next, looking at the region Voo < Vc, we see crox.p(ros) 
(circles) and crox.p(non-rcs) (triangles) plotted in green, 
with the corresponding analytically predicted result of 
eq. © (the two cross sections are expected to be roughly 
equal, from eq. 0) plotted as a solid line. crox_p(rcs) a-nd 
o'cx.p(non-res) appear to be roughly equal, with (Tcx.p(rcs) 
larger by ~ 10% for Voc <SC Vc, while the analytical pre- 
diction is roughly ~ 50% larger than the numerical re- 
sult. Plotted in red triangles is cr(;x_s(non-rcs)> along with 
the corresponding analytical result from eq. ((TJ plotted 
as a solid line. Clearly the agreement is superb. Plot- 
ted in red circles is the resonant cross section crox.s(ros): 
with the analytical result from eq. (|10|l plotted as a solid 
line. Although the statistics on this very small cross 
section are not as good as the other cross sections for 
Voc < fc, it still appears to agree very well with the 
analytical prediction. Finally, plotted in black circles is 
o'prcs(res)- Although from eq. Hll() above we expect it to 
be equal to (Texj3(res)i it is larger by nearly an order of 
magnitude. The discrepancy is due to the operational 
definition of resonance we have used — namely that the 
mean-squared separation of all bodies makes more than 
one minimum, with a successive minimum counted only 
if the intervening maximum between it and the previ- 
ous minimum is at least twice the value of both minima 
IklcMillan fcJlut Clearly this definition is too 

liberal, and classifies some preservation outcomes as res- 
onant when they should not be classified as such. In fact, 
requiring that the mean-squared separation make more 
than two minima for a resonance decreases crpi.(,s(-,.(,s) by 
about an order of magnitude, bringing it into agreement 
with the analytical prediction and making it nearly equal 
to crox_s(ros)- Requiring that intervening maxima be four 
times (instead of twice) the value of the minima on ci- 
ther side for a successive minimum to be counted has no 
discernible effect on the results. 

In the region Vc < v^c < I'orb resonance is, of course, 
forbidden. The cross section crox_p(non-ios) drops off 
quickly above Voo = Wc, obeying roughly 

^cx,p(non-rcs) ^ ^ 10^ f \ (16) 

for the case Mi = M2 = M, and q = m/M = 10"^. The 
cross sections for ionization and exchange of the incom- 
ing star for the host star (plotted in blue and red, re- 
spectively) are of the same order of magnitude, and only 
slightly different from the order of magnitude estimate 
in eq. IjlSfl (shown in the dashed line). In particular, for 
the bulk of the region, aion ~ 2CTox.s(non-rcs), and fTjon 
is ^ 20% larger than the order of magnitude estimate, 
yielding 

CTion « 1.27ra^r— , (17) 
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Fig. 4. — Statistics on the final planetary systems formed via non- 
resonant excliange in which the planet's original host star is ejected 
("ex_s (non-res)"), as a function oivac/vc, for the case Mi = M2 = 
M and q = 10~^. The fraction of final planetary systems that 
have hardened (semimajor axis decreased), /hardi shown in the 
solid line. The average change in the binding energy of the plane- 
tary system, {AEt)/Eh = {GMm/2a' - GMm/2a) / {GMm/2a) = 
a{l/a') — 1, is shown in the dashed line. The vertical dotted lines 
arc the critical velocity and the orbital speed. The horizontal dot- 
ted line represents the value 0.5 for /hard a-id for (AEi,) / Ei,. The 
hard-soft boundary lies where /hard crosses 0.5, and {AE^)/E^, 
crosses 0, which clearly and robustly occurs at Voo ~ 0.75dc. 
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Fig. 5. — The quantity (a/a')^^, the inversely- weighted average 
of the ratio of the planetary system's final semimajor axis a' to its 
initial value a, as a function of Vac/vc for the planetary systems 
formed via non-resonant exchange in which the planet's original 
host star is ejected ("ex_s (non-res)"), for the case Mi = M2 = M 
and q = 10""^. The vertical dotted lines show Vc and iiorbi while 
the horizontal dotted line shows the value 1. 



and 

O'ex.s(non-rcs) ~ O.bTTa — , (l»j 

V 

again, for the case Mi ~ AI2 — M, and q — m/M — 
10^"^. However, as discussed above in sectional we ex- 
pect the test-particle limit to be reached at <; « 10~^, so 
the results here should be valid for all q < 10^^. 

4. PROPERTIES OF EXCHANGE SYSTEMS AND THE 
HARD-SOFT BOUNDARY 

Having determined the character of the scattering in- 
teractions in each region of parameter space shown in 
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Fig. ^ we now consider the properties of the final plan- 
etary system for the outcomes which yield a planetary 
system — namely resonant preservation, and exchange of 
the host star for the incoming star. Since resonance is 
only allowed for Voo < Vc, for simplicity we consider only 
the outcome "ex_s (non-res)", which has a significant 
cross section for all ■^oo- Shown in Figure^jare the statis- 
tical properties of the final planetary systems in the out- 
come "ex_s (non-res)" as a function of Vao/yc, for the case 
Ml = M2 = M,q = 10"^ and e = 0. The fraction of final 
planetary systems that have hardened (semimajor axis 
decreased), /hard; is shown in the solid line. The average 
change in the binding energy of the planetary system, 
{AEb)/Eb = {GMm/2a' - GMm/2a) /{GMm/2a) = 
a{l/a') — 1, is shown in the dashed line. The verti- 
cal dotted lines are the critical velocity and the orbital 
speed. The horizontal dotted line represents the value 
0.5 for /hard and for {AEb)/Eb. For Voo ^ Vc, every 
outcome results in the planetary system hardening, with 
(AEb) / E), w 0.4. This implies that the semimajor axis 
has shrunk to about 70% of its initial value, as shown in 
Figure |5l which plots (a/a')^^ as a function of Voo/ Ve- 
in the region Vc < Voo < Worb, the planetary system pre- 
dominantly widens as a result of the encounter, with only 
« 30% of systems hardening. The systems don't widen 
by much, on average decreasing their energy by about 
20%. This corresponds to an increase in the semimajor 
axis of about 25%, as shown in Figure|31 This is rather in- 
teresting, since in this region exchange is nearly as likely 
as ionization, implying that a significant fraction of sys- 
tems that undergo a strong dynamical encounter survive 
it. And those that survive suffer a comparatively small 
change in their semimajor axis, with a significant frac- 
tion actually hardening. Finally, in the region Voo > Vorh, 
since the cross section for ionization is several orders of 
magnitude larger than that of exchange, the statistics 
are too poor to say anything conclusive. Although the 
curves in Figures ^ and |S1 appear to begin to rebound, 
we expect that they actually go smoothly to /hard = 0, 
{AEb)/Eb = -1, and (a/a')"^ ^ oo. 

As described in section El there has been consider- 
able variance in the literature on what velocity demar- 
cates the hard-soft boundary for systems with a wide 
mass ratio. A look at Figure ^ reveals that (at least for 
Ml = M2), the hard-soft boundary — which lies where 
/hard crosses 0.5 and {AEb)/Eb crosses — robustly oc- 
curs at Woo ~ 0.75wc. Of course, simply stating that the 
hard-soft boundary lies at Voo ~ 0.75vc belies the fact 
that for Vc < Woo < ^'orb the likelihood of a planetary 
system surviving a strong dynamical encounter in some 
form is rather significant. 

Based on the understanding garnered from our anal- 
ysis thus far, it is useful to consider the fate of a plan- 
etary system subject to strong dynamical encounters in 
a stellar environment of number density n and velocity 
dispersion v^- The average relative velocity between ob- 
jects in such a system is roughly \/2vcr, but since we are 
interested only in orders of magnitude here we will take 
Voo = VtT- Since we know all relevant cross sections, we 
can perform a simple ^^nuv'^ estimate for the lifetime of 
a planetary system as a function of a, 

r=^— , (19) 
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Fig. 6. — The cha racteristic hfetime r, plotted in the quantity 
TTi/v^ (see eq. 1221 '). of a planetary system with mass ratio q = 
10"'^, host star mass M = Mq, and eccentricity e = in a stellar 
environment of number density n and velocity dispersion Vac , as a 
function of a/ac- The critical semimajor axis ac is as defined in 
section |1] after eq. 1221 . The two vertical dotted lines show a = ac 
and a = a^j-^,. The quantity a^^-i, is the value of a for which the 
orbital speed in the planetary system is equal to v^. 



which we do now. A useful way to think about the 
evolution of a planetary system subject to encounters 
is the following. Taking the masses to be given, for a 
given value of a, we know Vc and Worb from eqs. and 
10). Given the value of Voo fixed by the velocity dis- 
persion of the stellar environment in which the plane- 
tary system finds itself, we thus have Voo/vc and know 
in which section of parameter space the planetary sys- 
tem lies in Figure 13 In the region Vac < Vc, by far 
the most likely outcome of a strong dynamical interac- 
tion is exchange in which the planet is ejected. Thus 
r for this region can be calculated using the cross sec- 
tion 0-cx.p = 0-ex.p(rcs) + O-ox.p(non-rcs) • In the regiou 

Voo > VoTb the dominant outcome is ionization, so one 
uses (Tion in eq. H19() . In between, in Vc < v^o < I'orb, the 
outcomes "ion" and "ex_s (non-res)" are almost equally 
likely. From eqs. and HHJ, roughly Xi = 2/3 of 
strong dynamical interactions result in ionization, while 
the remaining = 1/3 result in non- resonant exchange 
of the host star for the incoming star. The systems that 
undergo exchange suffer only a small perturbation to 
their orbital size, as evidenced by Figure|5| As a changes, 
Vc and Worb scale as a~^/^ with Vc/vorh fixed. Thus, in ef- 
fect, as a changes, the quantity v^c/vc scales as a^^^ with 
the distance between Vc and Uorb held fixed in log-space, 
and so the planetary system simply occupies a new posi- 
tion on the X-axis in Figure |31 A planetary system lying 
in regions Voo < Vc or v^c > Worb is destroyed as a re- 
sult of an interaction. But for a planetary system in the 
intermediate region, it is destroyed a fraction Xi of the 
time, while the other fraction Xc of the time it survives 
with Voo/vc increased by a mere ^ 10% on average (from 
Figure I^J. Eventually systems that survive via many 
exchange interactions will find themselves in the region 
Voo > Vorh and be destroyed quickly by ionization. Since 
the average change in a due to these exchanges is rather 
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small, let us, for the sake of simplifying the calculation 
of T, assume that they survive indefinitely. The average 
lifetime in this region of parameter space, including both 
ionization and exchange, is thus 

T^n + x,n + xln + xln + ■■■ = -^V ' (20) 

where 



^(f^ion "t~ f^cx_s(non— rcs))'^oo 

Since all cross sections of interest here scale as cr = 
(Joa/w^ where ctq is a function of the masses, we can 
be quite general and plot 

vl^ (To \acj GMim ' 

where the critical semimajor axis is Cc = GMim/ ^v^, as 
a function of a/cc in FigurclHl Note that the sharpness of 
the transitions at a = Oc and a = Oorb is due to the use of 
the analytical cross sections in calculating the lifetimes. 
From the numerical cross sections in Figure |31 it is evi- 
dent that, in reality, the transitions are a bit smoother, 
although the general structure is the same. Clearly the 
overall trend is for the lifetime to scale as a~-^, with the 
coefficient varying by only a factor of a few among the 
different regions of parameter space. Using only the fact 
that the hard-soft boundary lies at Voc = Vc, one might 
expect the lifetime to drop markedly at a = ac- In fact, it 
increases slightly there, and docs not drop until a = Oorb, 
due to the large cross section for planetary system pre- 
serving exchange in the intermediate regime. We now 
consider the implications of this result. 

5. IMPLICATIONS 

The main implication of our results is that the min- 
imum size of planetary systems that will likely be de- 
stroyed as a result of dynamical interactions in dense stel- 
lar environments is larger than has previously been as- 
sumed. This limit can be read almost directly from Fig- 
ure for a given velocity dispersion, stellar density, and 
timcscale. For example, in a globular cluster like 47 Tuc, 
which has a central stellar density of ~ 10^ pc~'^ and ve- 
locity dispersion of ^ lOkm/s, planetary systems wider 
than 0.05 AU will have been destroyed after ~ 10 Gyr. 
The second implication is that the statement that the lo- 
cation of the hard-soft boundary is at v^c = Wc, although 
correct, when taken on its own can be misleading. It is 
true that planetary systems on the low-velocity side of 
the hard-soft boundary that survive interactions tend to 
harden on average. However, planetary systems are over- 
whelmingly likely to be destroyed in this regime, as can 
be read directly from the cross sections in Figurc|31 Simi- 
larly, although planetary systems tend to soften, on aver- 
age, for iiqo > Vc, it is not until Voo > I'orb that ionization 
dominates the outcome. In other words, there is a wide 
region in parameter space (from Voo = Vc to v^o = Voih) 
for which the probability of exchange in which a plane- 
tary system results is comparable that of ionization, and 
those resulting planetary systems soften by only a small 
amount (as shown in Figures 01 and . 

Our results have a bearing on previous studies in the 
literature in which authors have adopted inappropriate 
prescriptions for outcomes of scattering interactions with 



planetary systems. As an example. iBonnell et al.l l|2001f) 
calculated the evolution of a population of planetary sys- 
tems in stellar environments of differing density, using 
a Monte-Carlo approach in which planetary systems un- 
dergo binary-single interactions with field stars. As a re- 
sult of each interaction, the planetary system is assumed 
to be ionized if v^o > Vc, and to harden if v^o < Vc- How- 
ever, as our results show, for Voo < Vc those planetary 
systems undergoing strong interactions are destroyed, for 
Vc < < "Worb they have an increased probability of 
surviving, and for Voo > Vorh they are predominantly 
ionized. The result of the discrepancy is that most of 
the systems destroyed on the high-velocity side of the 
hard-soft boundary [vao = Wc) in these simulations will 
actually survive as planetary systems of some kind. The 
growth of a spike in the relative frequency of planetary 
systems smaller than Oc predicted in low velocity dis- 
persion systems (such as open clusters), should not oc- 
cur, with the distribution of semimajor axes dropping 
smoothly from Oc to florb- There is the possibility for 
a slight hump in the relative frequency of planets with 
ttc < a < Qoih, based on our Figure |H1 but it would not 
be nearlv as pronoun ced as what is shown in Figure 5 of 
iBonneh et ani|200H) . 

On the observational side, we find that of the hot 
Jupiter planetary systems exp ected to be de t ected in 
47 Tuc in the transit survey of lGilliland elin ^QQ^— 
which was sensitive to periods up to « 4d — essentially 
all are expected to survive dynamical interactions in the 
high density core of the cluster on timescales comparable 
to the age of the cluster. This can be read directly from 
FigureEl with a » 0.05 AU. Thi s bolsters the conclusion 
drawn by iWeldrake et al.l l)200,^ that it is not dynamics 
that is responsible for the lack of hot Jupiters detected in 
47 Tuc , but rather may be metallicity-dependent effects 
(|Fische r & Valcnti 2005). The Wcldrak e et al., study 
sampled regions outside the core and was sensitive to 
orbital periods up to w 16 d, so the hot Jupiter binaries 
they looked for would also not have been destroyed by dy- 
namical interactions due to the lower density there. How- 
ever, it is possible that they would have been destroyed 
early in the lifetime of the cluster if it went through a 
high-density phase early on, the possibility of which is 
discussed in IBonnell et al.l l)2001l) . 

Looking ahead, we now speculate on future possibili- 
ties for observations of planets in dense stellar environ- 
ments. The globular cluster Terzan 5 is unique in that 
it has super-solar metallicity. Any argument against the 
formation of planets in globular clusters due to metal- 
licity effects is null in Terzan 5. If we adopt the fre- 
quencv of hot Jupiters in the so lar neighborhood of 0.8% 
ijButler et al.ll2000t lBasrill2004|) . and use the estimated 
central number den sity of lO^pc"'^ and core radius of 
0.5 pc ()Harrisiri99fiD . we expect ~ 100 hot Jupiters in 
orbit around main sequence stars in the core of Terzan 
5. According to Figure for = lOkm/s, all of these 
systems would survive for a Hubble time, since the defi- 
nition of "hot Jupiter" implies an orbit period less than 
~ 10 d. Terzan 5 lies nearly in the Galactic disk, close 
to the Galactic center, making optical observations of it 
practically impossible. However, it is readily observed 
in radio wavelengths, and already many millisecon d pul- 
sars (MSPs) have been discovered in it (,Ransom et alJ 
l2005j) . There is thus the possibility of detecting planets 
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in orbits around MSPs in this cluster. The most natural 
formation mechanism involves exchanging the MSP for 
the planet's host star in a dynamical interaction. From 
Figurc^we see that for the velocity dispersion of Terzan 
5 planetary systems with a ~ 0.01-10 AU lie in the in- 
termediate regime, and are thus most likely to possi- 
bly have exchanged in a MSP. The formation of such 
systems is a balancing act between destruction (for the 
wider systems) and a long interaction timescale (for the 
smaller systems), with a "sweet spot" at a ^ 0.1 AU. Al- 
though the expected number of MSP-planet systems is 
small (probably ~ 1), we can predict that such a system 
would have a ~ 0.1 AU. There is also the possibility that 
planets could form around puls ars, but the frequ ency of 
formation is rather less certain l)Sigurdssonlll992(l . 

6. PSR B1620-26 

The planetary system PSR B1620-26 in the globu- 
lar cluster M4 consists of a Jovian-mass object in an 
« 20 AU orbit about a MSP~WD binary The planetary 
system sits right at the edge of the core of M4 in projec- 
tion, which has a central velocity dispersion of ~ 4km/s. 
Looking at Figure ^ we see that the planetary system 
lies comfortably in the "intermediate" regime, making 
it very likely that it was formed via a dynamical scat- 
tering interaction in which the NS-containing binary ex- 
changed into the planetary system (as one unit) in place 
of the planet's original host star (scenario "2" below). 
With this motivation, we now perform a detailed numer- 
ical analysis of the possible formation scenarios of PSR 
B1620-26. 

First, it is useful to list the physical parameters of the 
system. The pulsar has a spin period of ~ 11ms and 
mass « 1.35 Af©, while the WD has mass 0.34±0.04 Mq. 
Hubble Space Telescope data has confirmed that the WD 
is young, with an age of ~ 0.5 Gyr. The PSR- WD bi- 
nary has eccentricity ~ 0.025. The planet has mass ~ 1- 
3Mj, where Mj is the mass of Jupiter, and orbits with 
eccentricity ^ 0.3 about the MSP-WD binary. The in- 
clination between the inner and outer orbits is ^ 55° j^g**. 
The currently favored dynamical formation scenarios are 
the following. In what we will call scenario "1" , a binary 
composed of a WD and a NS undergoes an exchange in- 
teraction with a binary consisting a main-sequence star 
near the turn-off mass and a Jovian planet. In the en- 
counter the WD is ejected from the system with the nor- 
mal star exchanging into orbit with PSR, leaving the 
planet in a wide orbit about the star-PSR binary. The 
star then undergoes mass transfer as it goes through the 
red giant phase, recyclin g the PSR to millisecond periods 
l|Sigurdsson et al.ll2'003l) . In scenario "2" , a binary com- 
posed of a NS and a near turn-off main-sequence star of 
mass 0.8 Mq undergoes an exchange interaction with a 
binary composed of a main-sequence star and a Jovian 
planet. In the encounter the NS-star binary acts as a sin- 
gle unit and exchanges into the planet orbit, ejecting the 
planet's original host star. The star then undergoes mass 
transfer as it goes through the r ed giant phase, r ecycling 
the PSR to milhsecond periods l|Ford et al.ll200nD . As we 
discuss below, our expectation that scenario "2" is the 
much more likely one proves to be correct. 

6.1. Quantifying the Outcome Likelihoods 
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Fig. 7. — Quantifying the likelihood of the possible formation 
scenarios of PSR B1620-26: The semi-normalized rate i?/[n7r(ao + 
ai)^] (in units of km/s) for various outcomes of a scattering inter- 
action between a NS-MSl binary and a MS2-P planetary system 
as a function of ao, the semi-major axis of the NS-MSl binary. 
Points are as defined in the legend. The rate for scenario "2" 
("[[NS MSI] P] MS2") clearly dominates over scenario "1" ("[[NS 
MS2] P] MSI") for a wide range in ao- 



As in section 13 we use the Fewbody numerical toolkit 
to perform scattering interactions and calculate cross sec- 
tions for different outcomes. The interaction rate per 
planetary system is simply 



R = 



O 



(23) 



where f{voo) is the velocity distribution function, n is the 
number density of NS binaries, and a is the cross section 
for the outcome of interest. Rates presented below were 
calculated by discretizing the integral over the distribu- 
tion function and summing over a range of velocities wide 
enough to guarantee convergence. In particular, we typ- 
ically integrate over a range of velocities from 1 km/s to 
15 km/s, with the upper l imit set by the es cape speed of 
15.25 km/s of the cluster (|Beer et al.ll2004|) . For the dis- 
tribution function we assume a simple Maxwellian with a 
velocity dispersion of Voo = V^Vm where = 3.88 km/s 
is the velocity dispersion in the core of M4 and the factor 
of is due to the fact that we are considering relative 
velocities between pairs of objects. Our naming conven- 
tion for the objects is analogous to what we used above, 
with "[NS MSI] [MS2 P]" denoting the initial configura- 
tion of the two binaries before the scattering interaction. 
"MSI" in fact represents a WD in scenario "1" and a 
main-sequence star in scenario "2" ; however, we label 
it "MSl" in both cases since it is only the mass of the 
object that matters in the scattering interactions (the 
radius of the object plays only a minor role when direct 
physical collisions are considered). We label the semi- 
major axis of the NS-MSl binary oq, and that of the 
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Fig. 8. — The normalized cross section, a/TT{ao + ai)^, for PSR 
B1620-26 formation scenarios "1" (stars) and "2" (circles) as a 
function of ai/ao, where ag is the initial semimajor axis of the 
neutron star binary and ai is that of the planet binary. The cross 
sections are calculated over a large grid in ag-ai space, hence their 
multi- valued nature as a function of ai/ao. All interactions were 
performed with Voo = 5km/s since the velocity-dependent rate 
peaks at this value. It is clear from the figure that as long as 
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Fig. 9. — The cumulative distribution of the center-of-mass ve- 
locity of the putative PSR B1620-26 triple after the exchange in- 
teraction in scenario "2", for the cases ao = 0.1 AU (solid lines) 
and ao = 1 AU (dotted lines) for several values of Vac from 
Ikm/s (left-most) to 15km/s (right-most). In both cases we set 
ai = 5AU. The triple system appears to suffer a recoil of ~ 3km/s, 
which is clearly not enough to eject it from the cluster (requiring 
fa 15km/s), but may be enough to displace it significantly away 
from the core. 



MS2-P binary ai. Since the NS binary must be signif- 
icantly tighter than the planet binary, it is expected to 
behave dynamically as a single unit; thus the cross sec- 
tions for different outcomes are expected to scale with ai 
(our numerical results confirm this). We therefore vary 
only ao, from 0.1 AU to 1 AU (a typical range of values), 
and keep ai fixed at 5 AU. We fix ei, the eccentricity of 
the planet binary, at 0, while we draw eo, the eccentricity 
of the NS binary from a thermal distribution truncated 
at large cq such that neither object in the binary exceeds 
its Roche lobe. Figure [7| shows the semi-normalized rate 
i?/[ri,7r(ao + ai)^] (in units of km/s) for the different out- 
comes listed in the legend as a function of oq. The most 
likely outcome is "[NS MSI] MS2 P", i.e., destruction of 
the planetary system with preservation of the NS binary. 
The next most likely outcome is scenario "2" ("[[NS MSI] 
P] MS2"), confirming our expectations from above. The 
next is any outcome involving direct physical collisions 
between system members, which we lump together into 
one category for the sake of succinctness. The next is 
scenario "1" ("[[NS MS2] P] MSI") which is larger than 
"[[NS MSI] MS2] P" only for oq > 0.3 AU. 

To verify that our results are not dependent on the 
choice of oi, we now consider the relative probabilities for 
scenarios "1" and "2" over a large grid in ap-oi space. 
To simplify the presentation of the results wc plot the 
normalized cross section for each outcome, (T/7r(ao+ai)^, 
as a function of ai/aa, as shown in Figure |H1 hence the 
multi- valued nature of the cross sections. For simplicity 



all interactions were performed with Voo = 5 km/s since 
the velocity-dependent rate peaks at this value. It is 
clear from the figure that as long as ai > Oq, which is 
quite reasonable on physical grounds, scenario "2" is the 
most probable. For the physically most reasonable value 
of ai/ag ~ 10, scenario "2" is more probable by more 
than an order of magnitude over scenario "1" . 

Having confirmed the preferred formation scenario, we 
now consider the orbital parameters of the hierarchical 
triple formed. Figure El shows the cumulative distribu- 
tion of the center-of-mass velocity of the triple after the 
exchange interaction, for the cases uq = 0.1 AU (solid 
lines) and ao = 1 AU (dotted lines) for several values 
of Voo from Ikm/s (left-most) to 15 km/s (right-most). 
The triple system appears to suffer a recoil of ~ 3 km/s, 
which is clearly not enough to eject it from the cluster 
(requiring « 15 km/s), but may be enough to eject it 
from the core. Figures EH and show the cumulative 
distributions of outer semimajor axis, Oout, and eccentric- 
ity, Cout, respectively, with the same convention for the 
line types as in the previous figure. The outer semimajor 
axis is clearly not particularly dependent on either ao or 
Woo, with a peak at ^ 10 AU. The outer eccentricity, on 
the other hand, while not particularly dependent on Voo, 
is strongly dependent on ao- For ao = 0.1 AU the eccen- 
tricity distribution is roughly thermal, with the cumula- 
tive distribution oc e^^^ . For ao = 1 AU the eccentricity 
distribution is flatter, with more systems at lower eccen- 
tricity. This result suggests that of the two, ao = 1 AU is 
the favored value for the NS binary semimajor axis, since 
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Fig. 10. — The cumulative distribution of the outer semimajor 
axis, aout, for the putative PSR B1620-26 triple system formed in 
scenario "2" . Line conventions and initial conditions are the same 
as in Figure 1^ 



Fig. 12. — The cumulative distribution of the inner semi-major 
axis, Oin, for the putative PSR B1620-26 triple system formed in 
scenario "2" . Line conventions and initial conditions are the same 
as in Figure 1^ 



0.4 



0.2 




0.2 0.4 0.6 O.E 

Outer Binary Eccentricity 

Fig. 11. — The cumulative distribution of the outer eccentricity, 
Eout, for the putative PSR B1620-26 triple system formed in sce- 
nario "2" . Line conventions and initial conditions are the same as 
in Figure 151 



we know from observations that the outer eccentricity is 
rathe r low, with eout ^ 0.2 ±0.1 (|Sigurdsson fc ThorsettI 
l2nnl . In Figm-cini we show the cumulative distribu- 
tion of the inner binary semimajor axis, with the 
same convention for the line types as in the previous 
figures. Confirming our earlier assumption, the inner bi- 



nary clearly behaves as a single dynamical unit, with the 
distributions for the final values strongly peaked around 
their initial values, and with very minimal dependence 
on Woo- 

6.2. Evolution of the Hierarchical Triple Due to Mass 

Transfer 

In scenario "2" , the hierarchical triple formed via the 
dynamical interaction is not the planetary system we see 
today in M4, but merely its progenitor. In this section 
we calculate the effects on the triple due to mass transfer 
in the inner binary as the main-sequence star ascends the 
red giant branch, and apply the results to the ensemble of 
triples formed in our scattering interactions. We assume 
slow isotropic mass loss, and tha t during this proce ss the 
eccentricities remain invariant ijRasio et al.lll99^ . For 
the sake of simplifying the calculation we also assume 
that the inner and outer eccentricities are zero. Since 
the outer orbit carries so little angular momentum (~ 
0.4% of the total angular momentum of the triple system) 
we neglect angular momentum coupling between the two 
orbits and treat them as two separate binaries during the 
mass transfer. 

The angular momentum of the inner orbit is given by 

Lin = AfNsA/MSl f , (24) 

where Min = A/ns + Mmsi, -^ns and Mmsi are the 
masses of the NS and the WD progenitor star, respec- 
tively, and Oin is the semimajor axis of the inner binary. 
Similarly, the angular momentum of the outer orbit is 
given by 



in 



MinA'/p 



Ga„ 



Mr 



1/2 



(25) 
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where Mt = A/in + A/p, Mp is the mass of the planet, 
and Cout is the semimajor axis of the outer orbit. Taking 
the logarithm of both sides and differentiating, we get 



NS 



MSI 



A/, 



NS 



MSI 



and 



SMu, SMv 



La 



Mi, 



^ 1 Sain 
2 flin 

1 faout 

2 ftout 



1 5A/i„ 

ISMt 

2 AfT 



(26) 



(27) 



for differential changes in the inner and outer orbits. 

In the mass transfer process, we assume that the mass 
lost from MSI gets channeled into an accretion disk 
around the NS, and any mass lost from the system is 
lost from the accretion disk with the characteri stic an- 
gular momentum of the NS (see e.g. lPodsiadlowski et al.l 
I2002|) . This imphes that 



6Lu 



' Ml " 



IM-n 



where fiin is the orbital velocity of NS, given by 



r^in — 



GAhr 



1/2 



(28) 



(29) 



Similarly, as mass is lost from the inner binary in the 
mass transfer process the angular momentum of the outer 
binary will change according to 



5L„ 



a. 



tSMr , 



(30) 



where flout is the orbital angular velocity of the inner 
binary in the outer orbit, and is given by 

n.u.^(^Y\ (31) 



•^out / 



We parameterize the mass loss by defining f3 as the 
fraction of mass lost from MSI that is accreted by the NS, 
so that (5A//ns = -PSMusi and (5A/i„ = (1 - i3)SMmsi- 
Note that with our conventions SMmsi < 0. Combin- 
ing equations H24(l through H31() . converting S^s to dif- 
ferentials and integrating, we find (jPodsiadlowski et al.l 
l2(ffll : 



^"MSl 
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MSI 
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for < /3 < 1 , (32) 
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X exp 



M, 
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A/ns 



A4 



for /3 = , (33) 



A/t 



r2Hi ^ for all P . 



(34) 



We can apply the expressions in eqs. (|32|) through 
(|34|) to the ensemble of triples formed via scenario "2" 
to determine their properties after mass transfer. Fig- 
ure 1131 shows the cumulative distribution of the inner 
binary semi-major axis after mass transfer. The solid 
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Fig. 13. — The cumulative distribution of the inner binary semi- 
major axis of the putative PSR B1620-26 triple system after mass 
transfer. The solid lines are for the case ag = 0.2 AU — in other 
words, the NS-MSl semimajor axis was 0.2 AU before the scat- 
tering interaction. The dotted lines are for the case ao = 0.3 AU. 
In both cases we set ai = 5AU and Voo = 5km/s. For each set 
of the lines the rightmost represents completely non-conservative 
mass transfer {/3 = 0), while the leftmost represents completely 
conservative mass transfer (/3 = 1). 
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Fig. 14. — The cumulative distribution of the outer semi-major 
axis of the putative PSR B1620-26 triple system after mass transfer 
for the cases oq = 0.2 AU (solid lines) and ao = 0.3 AU (dotted 
lines). Line conventions and initial conditions are the same as in 
Figure QSl 
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lines are for oq — 0.2 AU — in other words, the NS~MS1 
semimajor axis was 0.2 AU before the scattering inter- 
action. The dotted lines are for the case qq = 0.3 AU. 
In both cases we set oi = 5AU and Voo = 5km/s. For 
each set of the lines the rightmost represents completely 
non-conservative mass transfer (/3 = 0), while the left- 
most represents completely conservative mass transfer 
{f3 = 1). Figure FhI shows the cumulative distribution 
of the outer semi-major axis for the cases ag = 0.2 AU 
and flo ~ 0.3 AU. Line types are as in the previous figure. 
Comparing the distributions of the inner semimajor axis 
after the scattering interaction in Figure El with those 
after mass transfer in Figure El we see that the effect 
of mass transfer is to create a spread in the distribution 
of order the initial value. From observations we know 
the period of the inner orbit to be 191 d, corresponding 
to a s emimajor axis of « 0.8 AU ()Sigurdsson fc ThorsettI 
I2004|) . From this fi gure we see that we can quite reason- 
ably create such a system with uq « 0.3 AU for a wide 
range in the mass transfer parameter (3. From Figure [T^ 
we see that the distribution of the outer semimajor axis 
after mass transfer is much broader. Although the peak 
in the distribution appears to occur around 5-10 AU, it is 
quite plausible to form a system with an outer semimajor 
axis of 20 AU, as we observe for the system today. 

7. SUMMARY AND CONCLUSIONS 

We have studied in detail dynamical scattering interac- 
tions between star-planet systems and other single stars, 
tabulating analytical cross sections for all possible out- 
comes of the scattering encounter, and confirming them 
with numerical experiments (as shown in Figure |2J). In 
the scattering of planetary systems (or wide mass ratio 
binaries in general) there are two characteristic veloci- 
ties in the problem (as shown in Figure the critical 
velocity Vc, which is the value of Voo for which the total 
energy of the three-body system is zero; and Vorh, the or- 
bital speed in the binary. For a mass ratio q < 10""^ the 
two characteristic velocities differ by more than an order 
of magnitude, opening up a region in parameter space 
in which the nature of the scattering process is funda- 
mentally different from what occurs in comparable mass 
encounters. For Vc < v^o < ^^orb two outcomes dominate 
the cross with roughly equal weight: ionization, which 
destroys the planetary system; and exchange of the in- 
coming star for the planet's original host star. Thus in 
this region — which does not exist for comparable mass 
systems — there is a significant probability that a plan- 
etary system undergoing a strong dynamical encounter 
survives, albeit with a new host star. Since in the past 
workers in the field have typically taken as the ve- 
locity above which planetary systems are predominantly 
destroyed, this implies greater survivability of planetary 
systems in dense stellar environments than has been pre- 
dicted in some cases, as the characteristic lifetime plotted 
in Figure shows. 

We have also numerically determined the location of 
the hard-soft boundary for planetary systems, and find 
that it lies at v^o ~ O.TSwc. In other words, interac- 
tions with Voo < Vc in which the binary survives (e.g., 
via exchange) typically harden (a decreases), while for 
'^oo > Vc surviving binaries soften (a increases). How- 
ever, this simple statement belies the rich structure of 
the parameter space. In particular, for Vc < v^c < i^orb, 



although the surviving binary typically softens, it softens 
by only a small amount, increasing its semimajor axis by 
only ~ 25%, as shown in Figure |Sl Its final eccentricity, 
of course, is uncorrelated with the initial value, taking on 
a thermal distribution {dP/de = 2e). Only for Voo > Vovh 
does the binary soften by a large amount — that is, if it 
is not destroyed completely, which is the overwhelmingly 
dominant outcome for this region of parameter space. 

We have briefiy explored the implications of our re- 
sults. As mentioned above, the main implication is that 
the minimum size of a planetary system that will likely be 
destroyed as a result of dynamical interactions in dense 
stellar environments is larger than has been previously 
been assumed. For a stellar environment of a given ve- 
locity dispersion and density this limit can be read almost 
directly off of FigureEl Applying our results to the obser- 
vational campaigns searching for hot Jupiters in 47 Tuc 
via stell ar transits, we fi nd th at the planetary systems 
that the iGilliland et alJ (|200nD observations were sensi- 
tive to would have in fact survived to the prese nt day. 
This b olsters the conclusion drawn by iWeldrake et alJ 
()2005j) that it is not dynamics that is responsible for the 
lack of hot Jupiters detected in 47 Tuc, but rather may 
be metallicity-dependent effects. Looking toward the fu- 
ture, we predict a MSP-planet binary may likely be de- 
tected via pulsar timing in the super-solar mctallicity 
cluster Terzan 5, with a semi- major axis of ~ 0.1 AU. 

Recognizing that the planetary system PSR B1620-26 
in the globular cluster M4 lies in the "intermediate" re- 
gion of parameter space, we have performed a detailed 
numerical study of the possible scenarios in which it may 
have formed. By calculating cross sections and interac- 
tion rates, and treating the effects of mass transfer in the 
inner binary, we conclude that the scenario in which a 
NS-MS binary exchanges as a single dynamical unit into 
the planetary system for t he plan et's original host star 
(similar to the sce nario of iFord et al. 200U), is strongly 
favored over that of lSigurdsson et al... (,,,200311) , in which the 
exchange interaction involves ejecting the original NS bi- 
nary companion, with the planet's host star exchanging 
into it. This statement is quantified in Figure |H1 for a 
wide range of initial binary semimajor axes. 

Our general analysis (sections El through does not 
consider the effects of direct physical collisions between 
objects (e.g., MS stars and planets), which may play an 
important role when the ratio of stellar size to binary 
semimajor axis is > 0.01 (iFregeau et a,].ll200l . However, 
our code accurately treats physical collisions, and we 
have allowed for collisions in our analysis of PSR B1620- 
26 in sectionEl as shown in the cross sections of Figured 
We have not looked at the ejection speeds of the planets 
in encounters which destroy the planetary system. This 
topic has already been studied i n detail using jV-bod y 
methods llSmith BoTme11ll20f)TllH;H^v Sharall200^. 

Since our results assume point-mass particles, they 
are applicable to systems of any mass, as long as the 
mass ratio in the binary is relatively small < 0.1, and 
the incoming object and the heaviest object in the bi- 
nary are of comparable mass. We have not yet explored 
the application to scenarios involving, e.g., intermediate- 
mass black holes and super-massive black holes, stars 
and intermediate-mass black holes, and planetesimals 
and planets. Our results may also have limited appli- 
cation to wide mass ratio binaries containing more than 
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one light body, as long as the total mass in light objects 
is significantly less than the host mass, and the orbits of 
the light objects are not so strongly coupled that they 
transmit perturbations on timcscalcs comparable to the 
orbital period. 
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